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Abstract 

We perform an analytical and numerical study of the phase transitions in three- 
dimensional Z{N) lattice gauge theories at finite temperature for N > 4 exploiting 
equivalence of these models with a generalized version of the two-dimensional vector 
Potts models in the limit of vanishing spatial coupling. In this limit the Polyakov 
loops play the role of Z(N) spins. The effective couplings of these two-dimensional 
spin models are calculated explicitly. It is argued that the effective spin models 
have two phase transitions of BKT type. This is confirmed by large-scale Monte 
Carlo simulations. Using a cluster algorithm we locate the position of the critical 
points and study the critical behavior across both phase transitions in details. In 
particular, we determine various critical indices, compute the helicity modulus, the 
average action and the specific heat. A scaling formula for the critical points with 
N is proposed. 
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1 Introduction 



The phase structure of three-dimensional (3d) pure Z(N) lattice gauge theories (LGTs) 
has been the subject of an intensive study for more than three decades. It is well known by 
now that the zero-temperature models possess a single phase transition which disappears 
in the limit iV — > oo [lj. Thus, the U(l) LGT has a single confined phase in agreement 
with theoretical results [2]. The deconfinement phase transition at finite temperature 
is well understood and studied for iV = 2, 3. These models belong to the universality 
class of 2d Z(N) spin models and exhibit a second order phase transition in agreement 
with the Svetitsky-Yaffe conjecture |3j. Much less is known about the finite-temperature 
deconfinement transition when N > 4. The Svetitsky-Yaffe conjecture is known to connect 
critical properties of 3d Z(N) LGTs with the corresponding properties of 2d spin models, 
if they share the same global symmetry of the action. It is widely expected, and in many 
cases proved by either analytical or numerical methods, that some 2d Z(N > 4) spin 
models (like the vector Potts model) have two phase transitions of infinite order, known 
as the Berezinskii-Kosterlitz-Thouless (BKT) phase transitions. Then, according to the 
conjecture, the phase transitions in some 3c? Z(N > 4) gauge models at finite temperature 
could exhibit two phase transitions as well. Moreover, if the correlation length diverges 
when approaching the critical point, these transitions should be of the BKT type and 
belong to the universality class of the corresponding 2d Z(N) spin models. 

The BKT phase transition is known to take place in a variety of 2d systems: certain 
spin models, 2c? Coulomb gas, sine-Gordon model, Solid-on- Solid model, etc. The most 
elaborated case is the 2d XY model [H EJ [6]. There are several indications that this 
type of phase transition is not a rare phenomenon in gauge models at finite temperature 
- one can argue that in some 3d lattice gauge models the deconfinement phase transition 
is of BKT type as well. A well known example is the deconfinement phase transition 
in the U(l) LGT. Indeed, certain analytical [31 El E] as well as numerical results [HI [TU] 
unambiguously indicate the BKT nature of the phase transitions- 
Many details of the critical behavior of 2c? Z(N) spin models are well known - see the 
review in Ref. [11]. The Z(N) spin model in the Villain formulation has been studied 
analytically in Refs. [121 LT31 CLU 13 HH [17]. It was shown that the model has at least 
two phase transitions when N > 5. The intermediate phase is a massless phase with 
power-like decay of the correlation function. The critical index rj has been estimated 
both from the renormalization group (RG) approach of the Kosterlitz-Thouless type and 
from the weak coupling series for the susceptibility. It turns out that ^(fi^) = 1/4 at 
the transition point from the strong coupling (high-temperature) phase to the massless 

1 It should be noted, however, that the numerical results of 10 point to a critical index r\ larger than 
its XY value by almost a factor of 2 for N t — 8. Therefore, the question of the universality remains open 
for this model. 
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phase, i.e. the behavior is similar to that of the XY model. At the transition point /3c 
from the massless phase to the ordered low-temperature phase one has T}(0) = A/N 2 . 
A rigorous proof that the BKT phase transition does take place, and so that the massless 
phase exists, has been constructed in Ref. [IS] for both Villain and standard formulations 
of the vector Potts model. Universality properties of vector Potts models were studied 
via Monte Carlo simulations in Ref. [T9] for N = 6, 8, 12 and in Refs. [20| I2T1 [22| 123] 
for N = 5,7, 17. Results for the critical indices 77 and v agree well with the analytical 
predictions obtained from the Villain formulation of the model. 

We expect that 3d Z(N) gauge models at finite temperature exhibit the same critical 
properties as Z(N) spin models in two dimensions. In particular, gauge models with 
N > A may possess two phase transitions of the BKT type. On the basis of the Svetitsky- 
Yaffe conjecture the critical behavior of the gauge model in this case is governed by the 
2d Z(N) spin model. In particular, one expects the following values of critical indices: 
v = 1/2, r) — 1/4 at the first transition and v = 1/2, rj = A/N 2 at the second transition. 
To the best of our knowledge this scenario was not verified in the literature by either 
analytical or numerical means. The main goal of the present paper is to fill this gap and 
to study the nature of deconfining phase transitions in 3d Z(N) LGTs. 

The fact that the BKT transition has infinite order makes it hard to study its properties 
using analytical methods. In most of the cases studied one uses a renormalization group 
(RG) technique like in Ref. [12]. Unfortunately, there are no direct ways to generalize 
transformations of Ref. [12], leading to RG equations, to 3d Z(N) LGTs except for the 
limiting case N — > 00. To study the phase structure of these models we need numerical 
simulations. Here, however, another problem appears related to the very slow, logarithmic 
convergence to the thermodynamic limit in the vicinity of the BKT transition. It is 
thus necessary to use both large-scale simulations and combine them with the finite-size 
scaling methods. For a full finite-temperature gauge model this is an ambitious program, 
especially if one wants to study several values of N. We have therefore decided to utilize 
an approach developed by some of the present authors in Refs. [9j [10] and to divide 
the whole investigation into two steps. The finite-temperature model is formulated on 
anisotropic lattice with different spatial and temporal couplings. Here, following [9J, as a 
first step, we study the limit of vanishing spatial coupling. Since this approximation does 
not affect the global symmetry properties of the model, we hope that in this limit the 
model belongs to the same universality class as the full theory. It is well known that for 
both three- and four-dimensional gauge models at finite temperature this is indeed the 
case, at least for N = 2,3,4. Furthermore, in this limit the gauge model can be mapped 
onto a generalized Z(N) spin model in 2d. Thus, this investigation sheds some light on 
the details of the critical behavior of the general Potts model, i.e. beyond the vector 
Potts model, which is usually considered in the literature. In the present paper we study 
the phase transitions in models with N = 5, 7 in great details and in addition we locate 
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critical points of the models with N — 9, 13. Our computations are performed on lattices 
with temporal extent N t — 2, 4 and with spatial size in the range L G [128 — 2048]. 

This paper is organized as follows. In Section 2 we formulate our model and study 
some of its properties analytically. In particular, we establish the exact relation with a 
generalized 2d Z(N) spin model and discuss some of the RG predictions regarding the 
critical behavior. Also, we give simple analytical estimates for the critical couplings. In 
Section 3 we present the setup of Monte Carlo simulations, define the observables used in 
this work and present the numerical results of simulations for N = 5,7. In Section 4 we 
present results for the critical points in the models with TV = 9, 13 and discuss the scaling 
of the critical points with N. Our conclusions and perspectives are given in Section 5. 



2 Analytical considerations 

2.1 Relation of the 3d Z{N) LGT to a generalized 2d Z(N) vector 
model 

We work on a 3d lattice A = L 2 x N t with spatial extension L and temporal extension 
N t ; x = (xo,xi,X2), where xo G [0,N t — 1] and x±,X2 G [0, L — 1] denote the sites of the 
lattice and e n , n — 0, 1, 2, denotes a unit vector in the n-th direction. Periodic boundary 
conditions (BC) on gauge fields are imposed in all directions. The notations pt (p s ) stand 
for the temporal (spatial) plaquettes, l t (l s ) for the temporal (spatial) links. 
We introduce conventional plaquette angles s{p) as 

sip) = s n (x) + s m (x + e„) - s n (x + e m ) - s m (x) . (1) 

The 3d Z(N) gauge theory on an anisotropic lattice can generally be defined as 



Z(A;A,&;A0 = ]J ( 1 £ ] l[Q(s(p s )) l[Q(s(p t )) 

leA \ s(l)=ol Ps pt 



(2) 



The most general Z(A)-invariant Boltzmann weight with N — 1 different couplings is 



Q(s) = exp 



'JV-l 



k=i 



2nk 

cos s 

N 



(3) 



The Wilson action corresponds to the choice /3 P (1) = f3 p , /3 p (k) — 0, k — 2, N — 1. The 
U(l) gauge model is defined as the limit N — > oo of the above expressions. 
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To study the phase structure of 3d Z(N) LGTs in the strong coupling limit (f3 s = 0) 
one can map the gauge model to a generalized 2d spin Z(N) model with the action 



2 JV-1 



S 



x n=l k=l 



("■os | (s(x) - s(x + e n )) 



(4) 



The effective coupling constants ftk are derived from the coupling constant (3 t = /3 of the 
Z(N) LGT, using the following equation (the Wilson action is used for the gauge model): 



0k 



i f>(Q p ) cos 



(5) 



where 



Qk 



Bi 



N-l 

E 

p=0 
N-l 

J^exp 

p=0 



B ( 



f2npk\ 



(3 cos 



f27ipk\ 



(6) 
(7) 



These equations can be easily obtained in a few steps following similar computations for 
3D SU (N) model in the same limit f3 s = (see, for example, Ref. [23] and references 
therein) : 

• Fourier expansion of the original Boltzmann weight; 

• Integration over spatial gauge fields; this leads to an effective 2c? model for the 
Polyakov loops with the Boltzmann weight Q(s(x) — s(x + e n )) defined in ([6]); 

• Exponentiation and re-expansion in a new Fourier series. 

2.2 Renormalization group and critical behavior 

To gain some information on the critical behavior of the gauge model at (3 S = 0, one 
can perform the RG study of the theory. Let us consider the 2d model obtained after 
integration over spatial gauge fields, 



Z(A;ft AO 



n(^E) n 



. N 

xgA \ s(x)=0 



'N-l , 

(B k ) Nt cos ( 

.k=0 ^ 



2irk 
"AT 



(s(x) - s(x + e n )) 
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The coefficients B k = B k (f3) are given in ((?]) and can be represented as 



B k = tNr+kiP) • (9) 

r=—oo 

Here, I k (x) is the modified Bessel function. The spin variables s(x) can be associated 
with the Polyakov loops of the original model. 

The RG equations can be obtained only for the Villain formulation of the model. 
Replacing the Bessel function with its asymptotics and using the Poisson summation 
formula, one finds 




mftN) = iii^ > I ( 10 ) 



xGA \ s(x)=0. 



n 



m=— oo 



-^/3cff (j^( s ( x ) ~ s ( x + e n)) + 2irm^ 



This is nothing but the Villain formulation of 2d Z(N) vector model with /3 e ff defined as 

/3 eff = P/N t . (11) 

The RG equations for the model (fTTj) have been constructed in [12J. Their analysis has 
been performed by us in a recent work Therefore, all the conclusions as to the critical 
behavior remain valid in the present case. We shortly list the main results below: 

1. The critical RG trajectories in the planes ((3 e s,y) and ((3 e s,z), where y, z are the 
activities of the vortex configurations, coincide with those of the 2d Z(N) vector 
spin model; 

2. The critical index v has been computed numerically from the RG equations for 
a variety of N. It takes on the value 1/2 for all iV considered, in particular for 
iV = 5,7; 

3. The calculation of the two-point correlation function reveals that the index 77 is 
equal to 1/4 at the first transition point, while 77 = 4/iV 2 at the second critical 
point [12]. The same behavior is valid for the correlation of the Polyakov loops in 
our model; 

4. In |23j we have calculated the dependence of the critical points on N. These scaling 
formulae are expected to hold in the present case for any fixed N t and will be the 
subject of the discussion in Section 4. 
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2.3 Estimation of the critical points 



The location of the critical points in 3<i Z(N) LGTs for N > 4 is unknown. Therefore, 
before presenting numerical results, it is instructive to give some simple analytical predic- 
tions for the critical values of /3 crit for different values N t . Such a prediction could serve 
then as the starting point for the numerical search of the critical points. Following [9], 
such critical values can be easily estimated if one knows /3 cnt for N t = 1. Since the model 
with Nt — 1 coincides with the 2d Z(N) model, approximate critical points for other 
values of N t can be computed from the equality 



-DO 



-DO 



N t 



(12) 



where x on the right-hand side denotes the unknown critical point. Solving the last 
equation numerically one finds x. As for the critical points /3 crit at iV t = 1, we use 
our previous estimates from J2JJ. As it will be seen below, the predicted values are in 
a reasonable agreement with the numerical results. Table [1] summarizes approximate 
analytical predictions for the critical couplings at the two transitions, denoted by 8c and 
Pc , respectively; the last two columns show Monte Carlo values. 



Table 1: Values of BP and expected for N t = 1, 2, 4 in Z(N) with N = 5, 7. 



N 


N t 







Pc MC 


B (2) 

Pc MC 


5 


1 






1.051(1) 


1.105(1) 


5 


2 


1.8393 


1.9057 


1.87(1) 


1.940(7) 


5 


4 


2.7761 


2.8515 


2.813(3) 


2.898(4) 


7 


1 






1.111(1) 


1.88(8) 


7 


2 


1.9861 


3.1303 


2.031(7) 


3.366(7) 


7 


4 


3.2995 


4.9669 


3.406(8) 


5.158(7) 



3 Numerical results 

3.1 Setup of the Monte Carlo simulation 

To study the phase transitions we used the cluster algorithm described in [21]. The model 
is studied on a square L x L lattice A with periodic BC. Simulations were carried on for 
N t = 2, 4, but could easily be done also for other values of Nt, since the parameter 
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N t appears only in the definition of the couplings (J5J). As original action of the gauge 
model we used the conventional Wilson action. For each Monte Carlo run the typical 
number of generated configurations was 10 6 , the first 5 x 10 4 of them being discarded to 
ensure thermalization. Measurements were taken after 10 updatings and the uncertainty 
on primary observables was estimated by the jackknife method combined with binning. 
We considered the following observables: 

• complex magnetization Ml = \Mi\e % ^ ', 

M L = Y,exp(^s(x)) ; (13) 



population Sl, 

N /max i= o,7v_i l\ 

Sl ~ n^i{ 1? n) ' (14) 

where rii is number of s(x) equal to i; 

real part of the rotated magnetization Mr = \Ml \ cos(Nip) and normalized rotated 
magnetization = cos(Nifj); 

susceptibilities of Ml, Sl and Mr. xi M \ Xl \ Xl i 

= l 2 «- 2 > - O 2 ) ; (is) 



Binder cumulants U { T M) and b[ Mr) . 



jj{M) 



(\M L \<) 



3(|Mt| 2 } 2 ' 

B <„„> = <|Afa-WI') 
■\M R -{M R )\ 2 ) 



helicity modulus T, 
where 



T = (e)-L 2 (3(s 2 ) , (17) 



and the notation (ij) x means nearest-neighbors spins in the x-direction. 



S 



Figure 1: Scatter plot of the complex magnetization Ml at (3=1.6, 1.9 and 2.5 in Z(5) on 
a 16 2 x 2 lattice. 



3.2 Determination of the critical couplings 

A clear indication of the three-phase structure emerges from the inspection of the scatter 
plot of the complex magnetization Ml at different values of ft: as we move from low to high 
ft, we observe the transition from a disordered phase (uniform distribution around zero) 
through an intermediate phase (ring distribution) up to the ordered phase (N isolated 
spots), as Fig. [1] shows for the case of Z(5) on a 16 2 x 2 lattice. 

The first and most important numerical task is to determine the value of the two critical 
couplings in the thermodynamic limit, ft^ and ft^\ that separate the three phases. To 
this aim we have adopted several methods, which we list here: 

• Methods for the determination of ftc : 

(a) locate the position ftpc (L) of the peak of the susceptibility Xl ^ of the complex 
magnetization \Ml\ on lattices with various spatial size and find ft^ by a fit with 
the following scaling function, dictated by the essential scaling: 

/5(i) = /?(!) + d r , (18) 

P (\nL + B^ 

taking v equal to 1/2; 

(b) estimate the crossing point of the curves giving the behavior of the Binder 
cumulant U L versus ft on lattices with different spatial size L or, alternatively, 
search for the value of ft^ which optimizes the overlap of these curves when they 
are plotted against (ft — ft^){\nL) 1 ^ , with v fixed at 1/2; 

(c) consider the helicity modulus T near the phase transition and define ftpJ(L) as 
the value of ft such that rj(ft) = l/(2nftT) = 1/4 on the lattice with spatial size 
L [25], then find ftjp by a fit with the scaling function 
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Figure 2: Behavior of Ml and of its susceptibility versus f3 in Z(7) on lattices with N t = 2 
and L ranging from 128 to 2048. 



valid under the assumption that the phase transition belongs to the XY universality 
class. 

(2) 

• Methods for the determination of flc . 

(d) same as the method (a) using instead the susceptibility x L °f t ne population 
Sl', 

(e) same as the method (b) using instead simultaneously the Binder cumulant B^ 1 ^ 
and the order parameter m^. 

As an illustration of the methods (a) and (d), we show in Figs. [2] and [3] the behavior 
of Ml (and of its susceptibility) and that of Sl (and of its susceptibility) versus (3 in Z{7) 
on lattices with N t = 2 and L ranging from 128 to 2048. In Table [2] we summarize all the 
values of (3pJ(L) and j3$(L) found in this work for the application of methods (a) and 
(d) in Z(N) with iV = 5, 7 for N t = 2, 4. 

As an illustration of the method (b), we show in Fig. H]the behavior of versus 
(3 in Z(7) on lattices with N t = 2 and L ranging from 128 to 2048. Similarly, as an 
illustration of the method (e), we show in Figs. [5] and [6] the behavior of b[ Mr ^ and of 
versus (3 in Z(5) on lattices with N t = A and L ranging from 128 to 2048. 

As an illustration of the method (c), we show in Fig. [7] the behavior of the helicity 
modulus T versus (3 along with the line T = l/(2ir/3rj), r\ = 1/4, describing pseudocritical 
points, in Z(5i) and in Z(7) on lattices with N t = 2,4 and several values of the spatial 
extension L. For larger values of L only the points around the intersection were simulated. 

Finally, we report in Table |3] the determinations of the critical couplings (3^ and (3^ 
in Z(N) with iV=5 and 7 for N t =2 and 4, specifying the adopted method. 
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Figure 3: Behavior of Sl and of its susceptibility versus (3 in Z(7) on lattices with N t 
and L ranging from 128 to 2048. 
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Figure 4: Binder cumulant U L as function of f3 (left) and of (/3 — /3 c )(lnL) 1 / z/ (right 
Z(7) on lattices with N t = 2 and L ranging from 128 to 2048. 
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Table 2: Summary of all the determinations of /3 pc and (3p C in Z(N) on lattices with 
I? x N t . 



AT 
iV 


A7 


L 


R 1 

Ppc 


a 2 
Ppc 






128 


l. 7656(1) 


1 ^""70 f r* \ 

1.9773(5) 






192 
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1.9740(3) 






256 


1 7n 1 a/i \ 
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1.80124(7) 
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(Do 
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2.4913(9) 








32 


2.5928(6) 








64 


2.6538(5) 








128 


2.692(1) 


2.9376(8) 
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Figure 5: Binder cumulant B\ R function of (left) and of {(3 - /3 c )(ln L) 1 > v (right) 
in Z(5) on lattices with N t = 4 and L ranging from 128 to 2048. 
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Figure 6: as a function of (3 (left) and of ((3 — (3 c )(ln L) 1 ^ (right) in Z(5) on lattices 
with N t = 4 and L ranging from 128 to 2048. 
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Figure 7: Helicity modulus as a function of in Z(5) with N t = 2 (top left), in Z(5) with 
N t = 4 (top right), in Z(7) with N t = 2 (bottom left) and in Z(7) with N t = 4 (bottom 
right). 
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Table 3: Summary of the determinations of (3c and and (3c in Z(N) on lattices with 
size L 2 x N t . The fourth (sixth) column gives the method adopted to find (3^ {(3^)- 



N 


N 





method 


f3f ] 


method 






1.878(2) 


a 


1.91(6) 


d 


5 


2 


1.87(1) 


b 


1.940(7) 


e 






1.8801(7) 


c 










2.832(3) 


a 


2.88(3) 


d 


5 


4 


2.813(3) 


b 


2.898(4) 


e 






2.829(1) 


c 










2.070(2) 


a 


3.31(9) 


d 


7 


2 


2.031(7) 


b 


3.366(7) 


e 






2.069(1) 


c 










3.47(1) 


a 


5.0(2) 


d 


7 


4 


3.406(8) 


b 


5.158(7) 


e 






3.4782(8) 


c 







3.3 Determination of critical indices at the two transitions 

Once critical couplings have been estimated, we are able to extract some critical indices 
and check the hyperscaling relation. 

We start the discussion from the first transition. According to the standard finite-size 
scaling (FSS) theory, the equilibrium magnetization \Mj\ at criticality should obey the 
relation \Mj\ ~ L~^^ u , if the spatial extension L of the lattice is large enough H. Therefore, 
we fit data of \Ml\ at flip on all lattices with size L not smaller than a given L m ; n with 
the scaling law 

\M L \=AL-M"hi r L, (20) 

where a non-zero value for r takes into account the possibility of logarithmic correc- 
tions [261 EH- 

The FSS behavior of the susceptibility Xl ^ s gi ven by Xl ~ L 1 ^ , where 7/z/ = 2 — rj 
and 7] is the magnetic critical index. Therefore we fit data of Xl a ^ 011 au lattices 
with size L not smaller than a given L m i n according to the scaling law 

xi M) = AIT'" W L . (21) 

2 The symbol (5 here denotes a critical index and not, obviously, the coupling of the theory. In spite 
of this inconvenient notation, we are confident that no confusion will arise, since it will be always clear 
from the context which (3 is to be referred to. 
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As the value of the critical coupling (3^ we use the central value of the determination 
from the method (b) (see Table [3]). 

The results of the fits are summarized in Tables HI El El for the cases of Z(5) with 
N t =2, Z(5) with N t —4, Z(7) with N t =2 and Z(7) with N t =4, respectively. The reference 
value for the index 77 at this transition is 1/4, whereas the hyperscaling relation to be 
fulfilled is 7/1/ + 2/3 jv = d = 2. 

The procedure for the determination of the critical indices at the second transition is 
similar to the one for the first transition, with the difference that the fit with the scaling 
laws Eqs. (I2"0~j) and (12T|) is to be applied to data of the rotated magnetization, M R , and 
of its susceptibility, Xl* > respectively As the value of the critical coupling (3^ we use 
the one determined from the method (e) (see Table E])- 

The results of the fits are summarized in Tables El El HOI EH for the cases of Z(5) with 
N t =2, Z(5) with N t =4, Z(7) with N t =2 and Z(7) with N t =4, respectively. The reference 
value for the index 77 at this transition is 4jN 2 , i.e. 77 = 0.16 for N = 5 and 77 = 0.0816.. 
for N = 7, whereas the hyperscaling relation to be fulfilled is 7/1/ + 2(3 jv = d = 2. 

A general comment is that in many of the cases we investigated both d and rj at the 
two critical points slightly differ from the expected values, though these differences cancel 
to a large extent if we define rj as 2(3 jv. 

There is an independent method to determine the critical exponent 77, which does not 
rely on the prior knowledge of the critical coupling, but is based on the construction of 
a suitable universal quantity (2S1 EI]- The idea is to plot x^^' > L' n ~ 2 versus B{ MR) and 
to look for the value of 77 which optimizes the overlap of curves from different volumes. 
This method is illustrated in Fig. El for the case of Z(7) with N t = 4: for rj = 0.25 the 
overlap is optimal in the lower branch of the curves, corresponding to the region of the 
first transition, while per 77 = 0.0816 ~ 4/7 2 the overlap is optimal in the upper branch. 
Another option is to plot MrL 11 ' 2 versus m^, which leads to overlapping curves for 77 fixed 
at the value of the second phase transition, as illustrated in Fig. El for the case of Z(7) 
with N t = 2 and N t = 4. 
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Table 4: Critical indices (3/v and j/v for the first transition in Z(5) with N t = 2, deter- 
mined by the fits given in Eqs. (12~U|) and f[2"TT) on the complex magnetization Ml and its 
susceptibility XL a ^ A= = 1-87 for different choices of the minimum lattice size L min 
(an asterisk indicates a fixed parameter). The \ 2 of the two fits, given in the columns 
four and seven, is the reduced one. 



T 

-^min 


O 1 

P/v 


r fSlv 


2 

Xp/v 


7/1/ 


r-y/u 


2 

X-y/n 


a = 2p/v + j/u 


Tj = 2 — 7/^ 


16 


0.12040(7) 


0* 


4.43 


1.692(2) 


0* 


52.74 


1.933(2) 


0.308(2) 




0.1228(5) 


0.013(3) 


2.43 


1.97(1) 


-1.44(6) 


0.90 


2.22(1) 


0.03(1) 










1.668(2) 


0.125* 


61.75 


1.908(2) 


0.332(2) 


32 


0.12047(8) 


0* 


4.49 


1.712(2) 


0* 


22.49 


1.953(2) 


0.288(2) 




0.1239(7) 


0.019(4) 


2.12 


1.97(2) 


-1.4(1) 


0.99 


2.22(2) 


0.03(2) 










1.690(2) 


0.125* 


26.40 


1.931(2) 


0.310(2) 


64 


0.1206(1) 


0* 


4.30 


1.729(3) 


0* 


6.93 


1.970(3) 


0.271(3) 




0.126(1) 


0.034(7) 


1.51 


1.93(3) 


-1.2(2) 


0.73 


2.18(3) 


0.07(3) 










1.708(3) 


0.125* 


8.30 


1.949(3) 


0.292(3) 


128 


0.1208(1) 


0* 


2.43 


1.740(3) 


0* 


3.00 


1.982(3) 


0.260(3) 




0.126(2) 


0.03(1) 


1.72 


1.94(9) 


-1.3(5) 


0.84 


2.19(9) 


0.06(9) 










1.720(3) 


0.125* 


3.48 


1.962(3) 


0.280(3) 


192 


0.1210(1) 


0* 


2.09 


1.745(4) 


0* 


1.97 


1.987(4) 


0.255(4) 




0.125(3) 


0.03(2) 


2.06 


1.91(4) 


-1.1(2) 


0.97 


2.16(4) 


0.09(4) 










1.726(4) 


0.125* 


2.25 


1.968(4) 


0.274(4) 


256 


0.1211(2) 


0* 


2.19 


1.752(5) 


0* 


1.44 


1.994(5) 


0.248(5) 




0.124(4) 


0.02(3) 


2.56 


1.906(5) 


-1.02(3) 


1.19 


2.15(1) 


0.094(5) 










1.733(5) 


0.125* 


1.56 


1.975(5) 


0.267(5) 


384 


0.1212(2) 


0* 


2.44 


1.758(6) 


0* 


1.02 


2.000(6) 


0.242(6) 




0.1224(2) 


0.0079(9) 


3.27 


1.836(6) 


-0.53(3) 


1.37 


2.080(6) 


0.164(6) 










1.740(6) 


0.125* 


1.04 


1.982(6) 


0.260(6) 


512 


0.1215(3) 


0* 


2.18 


1.762(7) 


0* 


1.06 


2.005(8) 


0.238(7) 




0.1021(2) 


-0.135(1) 


1.21 


1.840(6) 


-0.54(3) 


1.88 


2.044(7) 


0.160(6) 










1.744(7) 


0.125* 


1.02 


1.987(8) 


0.256(7) 


768 


0.1209(4) 


0* 


1.71 


1.76(1) 


0* 


1.49 


2.00(1) 


0.24(1) 




0.1220(3) 


0.008(1) 


3.52 


1.834(8) 


-0.54(3) 


3.24 


2.078(8) 


0.166(8) 










1.74(1) 


0.125* 


1.47 


1.98(1) 


0.26(1) 


1024 


0.1209(6) 


0* 


3.41 


1.76(2) 


0* 


2.93 


2.00(2) 


0.24(2) 










1.74(2) 


0.125* 


2.89 


1.98(2) 


0.26(2) 



17 



Table 5: The same as Table I for Z(5) with N t = 4, determined at (3 ( c> = 2.813. 



-^min 






2 


7/1/ 




2 

Xj/ V 


a = 2p/v + "i/u 


r] = 2 — 7/1/ 


128 


0.1215(2) 


0* 


0.45 


1.730(5) 


0* 


1.23 


1.973(5) 


0.270(5) 




0.126(4) 


0.02(2) 


0.30 


1.808(6) 


-0.46(3) 


1.36 


2.06(1) 


0.192(6) 










1.709(5) 


0.125* 


1.32 


1.952(5) 


0.291(5) 


192 


0.1217(3) 


0* 


0.25 


1.734(6) 


0* 


1.20 


1.977(7) 


0.266(6) 




0.1228(3) 


0.007(1) 


0.32 


1.707(7) 


0.17(3) 


1.59 


1.952(7) 


0.293(7) 








0.25 


1.714(6) 


0.125* 


1.19 


1.957(7) 


0.286(6) 


256 


0.1218(3) 


0* 


0.30 


1.740(8) 


0* 


1.07 


1.983(8) 


0.260(8) 




0.1229(3) 


0.007(1) 


0.47 


1.62(3) 


0.7(2) 


0.88 


1.87(3) 


0.38(3) 










1.720(8) 


0.125* 


0.98 


1.963(8) 


0.280(8) 


384 


0.1217(5) 


0* 


0.44 


1.73(1) 


0* 


0.24 


1.97(1) 


0.27(1) 




0.1228(4) 


0.007(2) 


0.86 


1.785(9) 


-0.38(4) 


0.62 


2.03(1) 


0.215(9) 










1.71(1) 


0.125* 


0.22 


1.95(1) 


0.29(1) 


512 


0.1220(7) 


0* 


0.49 


1.72(2) 


0* 


0.05 


1.96(2) 


0.28(2) 










1.70(2) 


0.125* 


0.05 


1.94(2) 


0.30(2) 





1.6 




1.4 




1.2 




1 
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0.2 








r| = 0.25 
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Figure 8: Correlation between LP 2 and the Binder cumulant B^ Mr ^ in Z(7) with 
N t = 4 for i] = 0.25 (left) and for rj = 0.0816 (right) on lattices with L ranging from 128 
to 1024. 
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Table 6: The same as Table I for Z(7) with N t = 2, determined at 8 ( c > = 2.031. 



L ■ 


Blv 


T a 1 


yI, 




T I 


Y 2 , 


d = 28 lv + ilv 


71 = 2 — T / V 


1 98 
±ZO 


D 1 98fS9(Q N | 


n* 

u 


1 1 

-L -L . C/O 


1 7fi8(9"l 

L . 1 UOl Z 1 


n* 

u 


U.c/O 


9 09^(9^ 

Z.UZOl Z 1 


D 9^9(9") 

U.ZOZ I Z 1 




n i 41 (i 


u.u i o^y j 


U.c/O 


1 7fi^(9"l 

1 . i UO^Z J 


U.UOU ^O J 


1 1 ^ 


Z.U^O^O J 


fl 9^7(9") 

U.ZO i ^Z J 










1 747(9") 


fl 1 9^* 

U. ±ZO 


118 


9 nnzL(9'i 

z.uu^z J 


fl 9^(9"i 


1 Q9 


U.±Zt7Ul ± 1 


n* 

u 


7 7^ 
i . i o 


1 7fi7(^"l 

l./UI ^O J 


n* 

u 


1 DQ 


Z.UZO^O J 


U.iOO ^O J 




n izn(9"i 


n nod "\ 


n 87 

u.o < 


1 8Q(9"l 


n r(i "i 


n ^ 
u.oo 


9 1 7(^"l 


D 1 1 (9"l 

U. 11 ^Z J 










1 748(^ 


n 1 9^* 

U. ±ZO 


1 ^fi 


Z.UUU^O J 


fl 9 c >9(^'l 

U.ZOZ ^O J 




fl 1 9QzL(1 


n* 




1 771 (4") 


n* 

u 


fl 8^ 
u.oo 


9 D9Q(zn 
z.uzy I ^± J 


fl 99Q(4"l 

U.ZZc/ y-t ) 




n izl^c^ 


n in(9"i 

U. 1U ^Z J 


n q^. 


i.o / ^o j 


n r(9"i 

— U.U^Z J 


n 44 


9 1 fi(^ 

Z. lU^O J 


n 1 

U. lO ^O J 










1 752(4") 


125* 


1 00 


9 ni n(4"i 


D 948(4") 


384 


0.1297(2) 


0* 


4.96 


1.775(5) 


0* 


0.58 


2.034(5) 


0.225(5) 




0.1517(2) 


0.1480(8) 


0.27 


1.865(5) 


-0.61(2) 


0.55 


2.168(5) 


0.135(5) 










1.757(5) 


0.125* 


0.64 


2.016(5) 


0.243(5) 


512 


0.1303(2) 


0* 


2.43 


1.775(6) 


0* 


0.77 


2.035(7) 


0.225(6) 




0.1504(2) 


0.1389(9) 


0.38 


1.94(6) 


-1.2(4) 


0.32 


2.25(6) 


0.06(6) 










1.757(6) 


0.125* 


0.86 


2.017(7) 


0.243(6) 


768 


0.1309(4) 


0* 


0.76 


1.784(9) 


0* 


0.01 


2.05(1) 


0.216(9) 




0.1308(2) 


-0.001(1) 


1.54 


1.875(6) 


-0.65(3) 


0.00 


2.137(7) 


0.125(6) 










1.767(9) 


0.125* 


0.01 


2.03(1) 


0.233(9) 


1024 


0.1311(5) 


0* 


1.36 


1.79(1) 


0* 


0.00 


2.05(2) 


0.21(1) 










1.77(1) 


0.125* 


0.01 


2.03(2) 


0.23(1) 
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0.6 L = 1024 
L = 1536 
L = 2048 

0.4 



0.2 



r| = 0.0816 

Z(7), N, = 2 



0.2 



0.4 0.6 



0.8 



0.8 



0.6 



0.4 



L = 128 
L = 192 
L = 256 
L = 384 
L = 512 
L = 768 
L = 1024 



ti = 0.0816 



Z(7), N, = 4 



0.2 



0.4 



0.6 



0.8 



Figure 9: Correlation between MrLP' 2 and in Z(7) with N t = 2 (left) and in Z(7) 
with Nt = 4 (right) for 77 = 0.0816 on lattices with various values of L. 
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Table 7: The same as Table S] for Z(7) with N t = 4, determined at f3 ( c > = 3.406. 





P/v 


r P/» 


X/3/u 


j/u 


r-y/ u 




d = 2(3/v + -f/v 


rj = 2 — j/u 


128 


0.1294(1) 


0* 


4.64 


1.772(3) 


0* 


1.72 


2.031(4) 


0.228(3) 




0.141(2) 


0.07(1) 


0.54 


1.766(4) 


0.04(2) 


2.16 


2.049(9) 


0.234(4) 










1.751(3) 


0.125* 


1.76 


2.010(4) 


0.249(3) 


192 


0.1298(2) 


0* 


3.61 


1.771(5) 


0* 


2.12 


2.031(5) 


0.229(5) 




0.1456(2) 


0.0961(8) 


0.20 


1.85(8) 


-0.5(5) 


2.64 


2.14(8) 


0.15(8) 










1.751(5) 


0.125* 


2.20 


2.010(5) 


0.249(5) 


256 


0.1302(2) 


0* 


1.43 


1.778(6) 


0* 


1.84 


2.038(6) 


0.222(6) 




0.1432(2) 


0.0811(9) 


0.19 


1.74(3) 


0.2(2) 


2.58 


2.03(3) 


0.26(3) 










1.758(6) 


0.125* 


1.77 


2.018(6) 


0.242(6) 


384 


0.1307(3) 


0* 


0.37 


1.773(9) 


0* 


2.48 


2.03(1) 


0.227(9) 




0.1303(2) 


-0.003(1) 


0.78 


1.871(7) 


-0.63(3) 


5.32 


2.132(7) 


0.129(7) 










1.754(9) 


0.125* 


2.45 


2.01(1) 


0.246(9) 


512 


0.1310(5) 


0* 


0.08 


1.78(1) 


0* 


4.86 


2.04(1) 


0.22(1) 










1.76(1) 


0.125* 


4.77 


2.02(1) 


0.24(1) 



Table 8: Critical indices (3/v and j/u for the second transition in Z(5) with N t = 2, 
determined by the fits given in Eqs. ( )20l) and (12ip on the rotated magnetization Mr and 
its susceptibility Xl a ^ fic^ = 1-940 for different choices of the minimum lattice size 
L m i n (an asterisk indicates a fixed parameter). 



-^min 


P/v 




X 2 f3/u 








d = 2/3/v + -f/is 


rj = 2 — '-f/u 


128 


0.037(8) 
-0.27(3) 


0* 
-1.8(2) 


2.22 
1.07 


1.918(2) 
1.900(3) 


0* 
0.11(1) 


0.52 
0.53 


1.99(2) 
1.36(7) 


0.082(2) 
0.100(3) 


192 


0.03(1) 
-0.25(4) 


0* 
-1.6(2) 


2.10 
1.52 


1.918(3) 
1.821(3) 


0* 
0.59(1) 


0.64 
0.31 


1.97(2) 
1.33(8) 


0.082(3) 
0.179(3) 


256 


0.01(1) 
-0.20(5) 


0* 
-1.3(3) 


2.27 
2.38 


1.916(4) 
1.900(3) 


0* 
0.10(1) 


0.64 
0.82 


1.94(3) 
1.5(1) 


0.084(4) 
0.100(3) 


384 


-0.02(2) 
-0.01(2) 


0* 
0.03(6) 


1.77 
3.56 


1.913(6) 
1.87(3) 


0* 
0.3(2) 


0.62 
1.06 


1.88(5) 
1.85(6) 


0.087(6) 
0.13(3) 


512 


-0.05(3) 


0* 


1.37 


1.912(8) 


0* 


1.21 


1.81(7) 


0.088(8) 
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Table 9: The same as Table for Z(5) with N t = 4, determined at [3 ( c > = 2.898. 



Lm'm 


P/v 






j/u 






d = 2[3/v + y/v 


i] = 2 — j/u 


128 


0.17(1) 
-0.28(5) 


0* 
-2.6(3) 


1.52 
0.30 


1.850(3) 
1.847(4) 


0* 
0.02(2) 


0.07 
0.08 


2.20(3) 
1.3(1) 


0.150(3) 
0.153(4) 


192 


0.15(1) 
-0.1(1) 


0* 
-1.5(6) 


0.40 
0.25 


1.849(4) 
1.847(4) 


0* 
0.01(2) 


0.07 
0.09 


2.15(3) 
1.6(2) 


0.151(4) 
0.153(4) 


256 


0.15(2) 
-0.09(9) 


0* 
-1.5(5) 


0.51 
0.36 


1.849(5) 
1.846(5) 


0* 
0.01(2) 


0.07 
0.10 


2.14(4) 
1.7(2) 


0.151(5) 
0.154(5) 


384 


0.12(3) 
0.10(2) 


0* 

-0.15(9) 


0.24 
0.45 


1.847(8) 
1.845(6) 


0* 
0.01(2) 


0.07 
0.14 


2.09(7) 
2.04(5) 


0.153(8) 
0.155(6) 


512 


0.10(5) 


0* 


0.21 


1.85(1) 


0* 


0.11 


2.1(1) 


0.15(1) 



Table 10: The same as Table |8] for Z(7) with N t = 2, determined at AT = 3.366. 



-^min 


0/u 




Xp/v 


l/v 


r-y/ u 


Xy/y 


d = 2p/p + y/p 


1] = 2 — yjv 


128 


0.034(6) 
-0.27(5) 


0* 
-1.9(3) 


2.33 
0.34 


1.921(2) 
1.89(3) 


0* 
0.2(2) 


0.45 
0.27 


1.99(1) 
1.3(1) 


0.079(2) 
0.11(3) 


192 


0.018(7) 
-0.19(4) 


0* 
-1.3(3) 


0.72 
0.26 


1.919(2) 
1.903(2) 


0* 
0.10(1) 


0.23 
0.26 


1.96(2) 
1.53(9) 


0.081(2) 
0.097(2) 


256 


0.010(9) 
-0.12(7) 


0* 
-0.8(4) 


0.26 
0.18 


1.919(2) 
1.904(3) 


0* 
0.10(1) 


0.28 
0.28 


1.94(2) 
1.7(1) 


0.081(2) 
0.096(3) 


384 


0(1) 
-0.12(6) 


0* 
-0.9(4) 


11.54 
0.21 


1.919(3) 
1.904(3) 


0* 
0.10(1) 


0.35 
0.37 


2(2) 
1.7(1) 


0.081(3) 
0.096(3) 


512 


0.00(2) 
0.00(1) 


0* 
0.00(5) 


0.25 
0.38 


1.917(4) 
1.903(3) 


0* 
0.10(1) 


0.21 
0.29 


1.92(3) 
1.91(3) 


0.083(4) 
0.097(3) 


768 


-0.01(2) 
-0.01(1) 


0* 
0.02(6) 


0.11 
0.22 


1.916(6) 
1.903(4) 


0* 
0.10(2) 


0.29 
0.56 


1.90(5) 
1.89(3) 


0.084(6) 
0.097(4) 


1024 


-0.00(3) 


0* 


0.14 


1.918(8) 


0* 


0.53 


1.91(7) 


0.082(8) 
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Table 11: The same as Table El for Z(7) with N t = 4, determined at /3r = 5.158. 





P/u 






j/u 






d = 2p/v + j/v 


7] = 2 — 7/1/ 


128 


0.037(8) 
-0.27(3) 


0* 
-1.8(2) 


2.22 
1.07 


1.918(2) 
1.900(3) 


0* 

0.11(1) 


0.52 
0.53 


1.99(2) 
1.36(7) 


0.082(2) 
0.100(3) 


192 


0.03(1) 
-0.25(4) 


0* 
-1.6(2) 


2.10 
1.52 


1.918(3) 
1.821(3) 


0* 
0.59(1) 


0.64 
0.31 


1.97(2) 
1.33(8) 


0.082(3) 
0.179(3) 


256 


0.01(1) 
-0.20(5) 


0* 
-1.3(3) 


2.27 
2.38 


1.916(4) 
1.900(3) 


0* 
0.10(1) 


0.64 
0.82 


1.94(3) 
1.5(1) 


0.084(4) 
0.100(3) 


384 


-0.02(2) 
-0.01(2) 


0* 
0.03(6) 


1.77 
3.56 


1.913(6) 
1.87(3) 


0* 
0.3(2) 


0.62 
1.06 


1.88(5) 
1.85(6) 


0.087(6) 
0.13(3) 


512 


-0.05(3) 


0* 


1.37 


1.912(8) 


0* 


1.21 


1.81(7) 


0.088(8) 




Figure 10: Mean link energy (left) and specific heat (right) versus in Z(7) with Nt = 4 
on lattices with L ranging from 128 to 1024. 
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Table 12: Summary of the known values of the critical couplings (3c and (3 C in Z(N) 
with N t = 2, 4. 



N 


N t 






Reference 


5 


2 


1.87(1) 


1.940(7) 


this work 


7 


2 


2.031(7) 


3.366(7) 


this work 


9 


2 


2.04(3) 


5.38(4) 


this work 


13 


2 


2.02(1) 


10.815(8) 


this work 


oo 


2 




oo 




5 


4 


2.813(3) 


2.898(4) 


this work 


7 


4 


3.406(8) 


5.158(7) 


this work 


9 


4 


3.50(1) 


8.28(1) 


this work 


13 


4 


3.490(6) 


16.94(2) 


this work 


oo 


4 


3.42(1) 


oo 


m 



3.4 Other checks of the nature of the phase transitions 

To produce further evidence in favor of the fact that the phase transitions investigated so 
far are both of infinite order, we have calculated the mean link energy E and the specific 
heat C around the transitions in Z(5) and in Z(7) with N t =2 and N t = 4 (see Fig. [TOlfor 
the case of Z(7) with N t = 4, for example). In all cases the dependence of E and C on (3 
is continuous, so that first and second order transitions are ruled out. 

4 Behavior with N of the critical couplings 

The results of this work and those available in the literature allow us to make some con- 
siderations about the behavior with N of the critical couplings (3^ and (3^ ■ Examining 

(1 2) 

our data for [3 C , one concludes that, for a fixed N t , 

• (3c^ converges to the XY value very fast, like exp(— aN 2 ) , 

• (3c 2S> diverges like iV 2 . 

To better see the dependence on iV of the critical couplings, we have found also the 
critical coupling (3 C in Z(9) and Z(13) with N t = 2 and Nt = 4, using the method (b) 
for (3^ and the method (e) for [3^ (see Section 13731 for the details of these methods). In 
Table fT2l we summarize the present knowledge about the position of the critical points for 
3d Z(N) gauge models at (3 S = 0. 
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One should expect that the 3d Z(N) gauge models at (3 S = satisfy the scaling 
predicted by RG, probably up to O(N) corrections. One could try therefore to fit the 

fl 2) 

available Monte Carlo data for /3c ' with formulae predicted by RG and modified to 
account for such corrections. 

We find that the critical couplings for the first transition are well reproduced by the 
function 

/3«=A-5iV 2 exp (-^j , 

for suitable values of the parameters A, B and C, both for N t = 2 and N t = 4. 

The critical couplings for the second transition are well reproduced instead by the 
following functions: 

= AN 2 + BN + C 

and 

^ 2) = 1 ^-T^r+BN + C, 

1 - cos { w ) 

for suitable values of the parameters A, B and C, both for N t = 2 and N t = 4. The last 
formula has been suggested in Ref. [1] in the context of the zero-temperature theory. 

As an illustration of our fits, Fig. [11] shows the dependence of (3^ on iV (parameter 
B = in fitting formulas). One concludes from these plots that to distinguish between 
two scalings one should probably have data for smaller values of iV = 2, 3, 4 at (3 S = 0. 
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5 Conclusions 



In this paper we have studied the 3c? Z(N) gauge theory at the finite temperature in the 
strong coupling region (3 S = 0. This study was based on the exact relation of this model to 
a generalized 2d Z(N) spin model. In Section 2 we established the exact relation between 
couplings of these two models, described qualitatively some of the RG predictions for 
effective model and gave analytical estimations of the critical couplings. 

The main, numerical part of the work has been devoted to the localization of the 
critical couplings and to the computation of the critical indices: 

• We have determined numerically the two critical couplings of Z(N = 5, 7, 9, 13) 
LGTs and given estimates of the critical indices rj at both transitions. For the first 
time we have a clear indication that for all N > 5 the scenario of three phases is 
realized: a disordered phase at small (3 t , a massless or BKT one at intermediate 
values of (3 t and an ordered phase, occurring at larger and larger values of (3 t as N 
increases. This matches perfectly with the N — > oo limit, i.e. the finite-temperature 
3d U(l) LGT (at (3 S = 0), where the ordered phase is absent; 

• We have found that the values of the critical index 77 at the two transitions are 
compatible with the theoretical expectations; in order to reproduce the expected 
value of 77, we had to take into account logarithmic corrections with the exponent r 
fixed to 0.125; 

• The index v also appears to be compatible with the value 1/2, in agreement with 
RG predictions. 

Results listed above present further evidence that finite-temperature 3c? Z(N) LGTs 
for N > 4 undergo two phase transitions of the BKT type. 

Moreover, this model belongs to the universality class of the 2d Z(N) spin model, at 
least in the strong coupling limit (3 S — 0. 

Considering the determinations of the critical couplings as a function of N, we have 
conjectured the approximate scaling for (3c ' (N). 

Finally, the study performed here allows to improve our knowledge of the phase di- 
agrams of the generalized 2c? Z(N) spin models. As an example, we plot in Fig. [T21 the 
general phase diagram for N = 5 in the (ti, t 2 )-plane, where £, = (Bi/B ) Nt , and Bi are 
defined in (|7j). Here, the line AB is self-dual line, SPM corresponds to the standard Potts 
model, VPM to the vector model. The SPM undergoes a first order phase transition with 
a critical point occurring on the self-dual line. The line of the first order phase transition 
terminates at the point 2. Its approximate position was computed in Ref. |29j. Shown 
are also the locations of the critical points for the VPM (N t = 1) and for Z(5) LGT with 
N t = 2, 4. The parametric curves for different N t lie very close to each other, so we cannot 
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Figure 12: Phase structure of the general Z(5) spin model (see text for explanation). 

trace a sufficiently big part of the curve while changing N t . On the other hand it shows 
that already the model with N t = 4 presents a very good approximation to the finite- 
temperature limit. Indeed, the parametric curve for Nt = 8 is almost indistinguishable 
from the curve with N t = 4. 

6 Acknowledgments 

The work of O.B. was supported by the Program of Fundamental Research of the De- 
partment of Physics and Astronomy of NAS, Ukraine. The work of G.C. and M.G. was 
supported in part by the European Union under ITN STRONGnet (grant PITN-GA- 
2009-238353). G.C. thanks Matteo Giordano for useful discussions. 

References 

[1] G. Bhanot and M. Creutz, Phys. Rev. D 21 (1980) 2892. 

[2] A. Polyakov, Nucl. Phys. B 120 (1977) 429; T. Banks, J. Kogut, R. Myerson, Nucl. 
Phys. B 121 (1977) 493; M. Gopfert, G. Mack, Commun. Math. Phys. 81 (1981) 97. 

[3] B. Svetitsky, L. Yaffe, Nucl. Phys. B 210 (1982) 423. 



26 



[4] V. Berezinskii, Sov. Phys. JETP 32 (1971) 493. 
[5] J. Kosterlitz, D. Thouless, J. Phys. C 6 (1973) 1181. 
[6] J. Kosterlitz, J. Phys. C 7 (1974) 1046. 
[7] N. Parga, Phys. Lett. B 107 (1981) 442. 

[8] O. Borisenko, Critical behaviour of 3D U(l) LGT at finite temperature, PoS LAT 2007 
(2007) 170. 

[9] O. Borisenko, M. Gravina, A. Papa, J. Stat. Mech. 2008 (2008) P08009. 

[10] O. Borisenko, R. Fiore, M. Gravina, A. Papa, J. Stat. Mech. 2010 (2010) P04015. 

[11] F.Y. Wu, Rev. Mod. Phys. 54 (1982) 235. 

[12] S. Elitzur, R.B. Pearson, J. Shigemitsu, Phys. Rev. D 19 (1979) 3698. 
[13] M. B. Einhorn, R. Savit, and E. Rabinovici, Nucl. Phys. B 170 (1980) 16. 
[14] C.J. Hamer, J.B. Kogut, Phys. Rev. B 22 (1980) 3378. 
[15] B. Nienhuis, J. Statist. Phys. 34 (1984) 731. 
[16] L.P Kadanoff, J. Phys. A 11 (1978) 1399. 
[17] J.L. Cardy, J. Phys. A 13 (1980) 1507. 

[18] J. Frohlich, T. Spencer, Commun. Math. Phys. 81 (1981) 527. 
[19] Y. Tomita, Y. Okabe, Phys. Rev. B 65 (2002) 184405. 

[20] O. Borisenko, G. Cortese, R. Fiore, M. Gravina, A. Papa, Critical properties of the 
two-dimensional Z(5) vector model, PoS LATTICE 2010 (2010) 274 |arXiv: 1101 .05121 
[hep-lat]]. 

[21] O. Borisenko, G. Cortese, R. Fiore, M. Gravina and A. Papa, Phys. Rev. E 83 (2011) 
041120. 

[22] O. Borisenko, G. Cortese, R. Fiore, M. Gravina, A. Papa, Critical properties of 2D 
Z(N) models for N > 4, PoS LATTICE 2011 (2011) 304 |arXiv:1110.6"385l [hep-lat]]. 

[23] O. Borisenko, V. Chelnokov, G. Cortese, R. Fiore, M. Gravina, A. Papa, Phys. Rev. 
E 85 (2012) 021114. 



27 



[24] M. Billo, M. Caselle, A. D'Adda, Int. J. Mod. Phys. A 12 (1997) 5753. 

[25] A. Yamagata and I. Ono, J. Phys. A24 (1991) 265. 

[26] R. Kenna and A. C. Irving, Nucl. Phys. B 485 (1997) 583. 

[27] M. Hasenbusch, J. Phys. A 38 (2005) 5869. 

[28] D. Loison, J. Phys.: Condens. Matter 11 (1999) L401. 

[29] E. Domany, D. Mukamel, A. Schwimmer, J. Phys. A 13 (1980) L311. 



28 



